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RESONANCE ET FORCES DE GRAVITATION. | 


. Depuis Enr a signalé, voici bientôt 50 ans, la possibilité d’utiliser les 
effets de résonance pour amplifier les oscillations de la bas de torsion sous 
l’action des forces de gravitation, il ne me semble pas qu’on en ait tiré un parti 
suffisant, soit pour réaliser d’une maniére simple la mise en évidence de ces 
forces devant un auditoire, soit pour augmenter la précision obtenue jusqu'ici 
dans la mesure de la constante de Newton. 

„Soient deux systèmes matériels rigides A et B mobiles autour d’un même axe 
vertical et possédant l’un et l’autre un plan de symétrie passant par cet axe. 
Le système A est suspendu, suivant l’axe, à un fil de torsion. Nous prendrons 
pour origine des azimuts la position du plan de symétrie de À lorsque ce système | 
est en équilibre en l’absence de B. Pour toute autre orientation de A, l’azimut . 
du plan + symétrie sera désigné par a; soit 9 l’azimut du plan de symétrie de B. 
Par raison de symétrie, le moment C par rapport à l’axe des forces de gravitation 


“exercées par B sur A ne dépend que de l’écart angulaire 0 — a et en est une fonc- 


tion périodique impaire de période 27 si A et B ne possèdent pas d’autres plans 
de symétrie passant par l’axe. Le développement de cette fonction en série de 
Fourier est donc de la forme : 


(1) C=C, sin (f — a) +C, sin 26 — a) + ..- +C, sinn(@—a)-+... 


Il est facile, comme nous le verrons sur un exemple, de disposer les masses 
de A et de B de manière que les coefficients G,, G,... forment une série rapi- 
dement décroissante. 

Nous supposerons qu'on donne au système B un mouvement connu, par, 
exemple une rótation uniforme de vitesse angulaire w, autour de l’axe commun. 
Le mouvement de A est alors régi par l'équation : 


| da da WER . Er 
(2) 175 HA q +Ka= 20. sin n (wt — a) 
en supposant qu’à l’origine du temps le plan de symétrie de B coïncide avec 
Pazimut origine. I est le moment d'inertie de A par rapport à l’axe commun, 
à le coefficient qui determine l'amortissement des oscillations de À, et K le coeff- 
cient de torsion du fil de suspension, modifié, s’il y a lieu, par les autres actions 
directrices auxquelles peut être soumis le système A, Comme terme de compa- 
raison pour la déviation a, nous prendrons la quantité : 

C, 

(3) GE 


que nous appellerons la déviation statique et qui coïncide très sensiblement, pour 
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les systèmes à C, rapidement décroissants, avec la déviation statique de A lors- - 
que l'écart angulaire 0 — a prend la valeur, voisine de F’ pour laquelle le 


moment C est maximum. La petitesse des actions de gravitation fait que l'angle - 
ay sera généralement petit, de l’ordre d'un degré au maximum. 

Si le but poursuivi, celui de démonstration expérimentale, par exemple, 
n'exige pas une grande précision, nous pourrons nous contenter d'une pre- 
mière approximation dans laquelle nous négligerons, au second membre de (2), 
la déviation a devant wt. 

L’équation (2) devient alors Pure et la solution a est la somme de termes a, 
qui correspondent respectivement aux différents termes de la somme qui 
figure au second membre. 

Considérons tout d’abord la solution de régime permanent; elle aura pour 


expression, en désignant par Q la pulsation vr des oscillations propres du 


systeme A : 
(4) | ~ a = ŽA, sin (nwt + ®,) 


avec : 
C, 


ae on mp) 
2,202 
qye a+ PEE 
T 


nwQd 
tg Pn = ro 


Aie 


ô étant le décrément logarithmique : 


RY 
=i 
des oxcillatioós propres du système À. 

Si l’on donne à la rotation uniforme de B une vitesse angulaire w précisément 
égale à la pulsation propre O des oscillations de A, on obtient un effet de réso- 
nance dans lequel le premier terme de (4) l'emporte de beaucoup sur les autres 
en amplitude et prend la forme simple : 


ai= = cos Ot 
3 pouvant être rendu très petit, on voit que l'amplification 7 de l'amplitude 


de a, par rapport à la déviation statique, peut être rendue ès grande. 

Si Pon se fixe une amplification donnée.A pour l'amplitude des oscillations 
de régime permanent par rapport à la déviation a, le décrément logarithmiqne 
se trouve fixé par : 


dach‘ 
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ainsi que le temps nécessaire pour l'établissement de cette amplitude á une 


approximation donnée, au centième près par exemple. On obtient facilemént 
pour ce temps : 


A 
lo 


Au point de vue de la rapidité des mesures ou de la démonstration expéri- - 
mentale, il est avantageux de s'adresser au régime variable de préférence au 
régime permanent, en donnant au système À un amortissement aussi faible que 
possible. Nous supposerons qu’à l’origine du temps les angles a et 6 sont 
nuls et que l'équipage A est immobile dans sa position d’équilibre, la valeur 


TES XI 5 . . A . 
initiale de se trouvant ainsi nulle. Le terme d'amortissement dans le premier 


membre de l’équation (2) pouvant être considéré comme négligeable au début 
du régime variable, on obtient pour la solution «, qui correspond en première 
approximation au premier terme du second membre, et à la résonance : 


(6) | a =? [sin Qt — Ot cos Qi). 


Les solutions a,... %, qui correspondent aux autres termes de la somme 
qui figure au second membre de l'équation (2) sont toutes de la forme : 


x 


C, . . ~ 
Os Re) [n sin Ot — sin ni]. 


Il est remarquable que, pour les valeurs de t égales aux multiples entiers A 
de la période = des oscillations propres de A ou de la rotation de B, les a, 


s’annulent tous, à l’exception de a, qui donne : 
(7) a, = kray. 


L’amplification A égale à kr que représente cette lem: élongation négative 
par rapport à la déviation statique a, est ainsi obtenue au bout d'un temps : 


La comparaison avec l'équation (5) montre l'avantage que présente, sur 
l'observation du régime permanent, l'observation des élongations négatives en 
régime variable de. résonance. Il faut y ajouter l'avantage que présente rél- . 
mination des termes autre que celui en C, dans le second membre de 1 équa- 
tion (2). i i 

-En utilisant la disposition classique de M. Vernon Boys avec une suspension 


22. 


Y 
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robuste et une durée d'oscillation d'une minute, on obtient facilement une 
déviation statique d'un demi-degré; la formule (7) montre qu'après quelques 
minutes de rotation de B, l’élongation négative a prend une valeur aussi 
grande qu’on peut le désirer pour une observation facile à distance. 

Si Pon a en vue la détermination de la constante G de Newton, l'observation 
des élongations négatives en régime variable de résonance me semble permettre, > 
dans des conditions plus commodes, d'augmenter la précision obtenue jus- | 
Ce par l'observation de la déviation statique dans des balances de torsion 

e trés longue durée d’oscillation. i . 

La mesure de G se raméne á la détermination, en valeur absolue, du moment 
C,. Celle-ci.peut s’obtenir par une méthode de zéro, en compensant C, par 
une action électromagnétique ou électrodynamique alternative, Pannulation 
de l’élongation négative, qui dépend seulement du terme fondamental dans le 
couple PR C exercé sur A, servant de critérium sensible pour cette 
compensation. On peut aussi procéder directement par observation ou enre- 
gistrement du régime variable tel qu'il est représenté, en première approxima- | 
tion, par l'équation (6), avec des corrections qu'il est facile de déduire de la Ñ 
méthode suivante pour l'intégration de l'équation (2) par approximations 
successives. En divisant Jes deux membres de cette équation par K nous la 
mettrons sous la forme : $ 


1 de de à e i + 
(9) or tra tem aulin@t—a)+e, sin 2wt— a) + ...] 


_en posant : 


O | | 
y Le = ü, | 

Les nombres ¢,, ..., €, sont petits et décroissent rapidement quand n 
augmente. Comme ils interviendront seulement dans des termes correctifs, on y 


les déterminera avec une précision suffisante en appliquant la méthode des 
élongations négatives aux régimes variables de résonance sur les harmoniques = 
successifs de la pulsation fondamentale w. En faisant tourner B avec la vitesse | 


angulaire a, c'est le ni°m° terme de la somme qui figure au second membre, 


de l'équation (2) qui se trouvera en résonance avec Poscillation du système A, 
et les élongations négatives en régime variable ne dépendront que de C,. 
L’angle ap, mesuré en radians, étant très petit, nous obtiendrons la solution a _ 
de l'équation (9) sous forme d'une série très rapidement convergente en la 
développant suivant les puissances croissantes de a, sous la forme : 


aa AA Ea faft) + + + 


Les équations auxquelles satisfont les fonctions f s’obtiennent en substituant 
à a dans (9) le développement (10) et en égalant les coeflicients des mêmes 


i | 
7 


v 


4 
d 
3 
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puissances de a, dans les deux membres. On obtient ainsi, en supposant la réso- 
nance sur le terme fondamental : 


i@f i yd : 
aie te = sin Ot+c, sin20t+... +e, sin nQt, 


Chaque second membre F, s'exprime en fonction des f antérieurement 
connus, depuis f, jusqu’à f,_;. Chaque fonction f est déterminée par la condi- 
tion de s’annuler à l’origine du temps ainsi que sa dérivée par rapport à t. 
Les intégrations des équations (11) se font finalement par des calculs numé- 
riques classiques qu'il n'est d’ailleurs pas nécessaire de pousser bien loin pour 
obtenir une approximation largement suflisante. 

On obtient ainsi l'équation dont (7) est la première approximation et qui 
donne les élongations successives en régime variable. De l'observation de ces 
élongations on déduira a, et par conséquent Cı. y 


Pour obtenir G, il faut connaître la valeur de C, pour une combinaison de 
deux systèmes A et B de structure assez simple pour que les actions de gravi- 
tation entre eux soient facilement calculables. Br 


On fera, par exemple, en sorte que l’équipage oscillant puisse se composer 
d’un système de sphères homogènes A et d'un support A’ ou de ce support 
seul, et que la partie tournante puisse de la même manière se composer d'un 
système de sphères homogènes B et d'un support B’ ou de ce support seul, 
Le caractère additif des actions de gravitation a pour conséquence la relation 
évidente : ` 


C,(A,B) = C,(A + A’,B + B') — C (A + A',B' | 
u) (A,B) = C(A + A,B +B) à Aa 


On déterminera expérimentalement, comme il a été dit, les quatre termes 
du second membre, ce qui donnera par (12) la valeur expérimentale de C,(A, B). 


Le méme caractére additif des actions de gravitation permet de calculer la 
valeur théorique de C,(A, B) en associant successivement chacune des sphéres 

ui composent À avec chacune de celles qui composent B. Soit, pour l’un 
iy ces couples, a la distance à Paxe du centre de la sphère A, b la distance à 
laxe du centre de la sphère B et h la distance des projections de ces centres 
sur laxe, A Pazimut du centre de la sphère Á par rapport au plan de symétrie 
du système oscillant, m la masse de cette sphère, B Pazimut du centre de la 


sphère B par rapport au plan de symétrie du système nis M la masse de 
a ö 


cette sphère. On obtient facilement en posant p = 
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; C GMm cos(A—B) 2% cos 6 dû 
13) es pora), VIE p cos 6 


2 ELE ja ale TFP Eg] 
D td Vit? V1- p E(k?) 


F(k?) et E (k?) sont les intégrales elliptiques complètes de Legendre pour le 
module défini par : = 
aa RE 
| £ itp 
La sommation qui figure dans (13) doit être étendue à toutes les combinai- 


sons d'une sphère de A avec une sphère de B. 


On obtiendrait de la même manière pour le coefficient C, du terme général 
de la série (1). 


c y GMmn 12%” cos nb dû AB 
(a) rl rr X 008 (AB). 


L'intégrale qui figure dans (14) se calcule également à partir des fonctions 
de Legendre.  : 


Quand les p sont petits, c’est-à-dire quand le système A est très ramassé autour 


de laxe de suspension, ce qui est d’ailleurs favorable à la sensibilité, le rapport : 


e, de €, à C, est de l’ordre de : 


(PL 3.5.7... (2n—1) 
a= (£) 246... (2n—2 


C'est-à-dire que les C, décroissent sensiblement, quand n augmente, comme 
les termes d'une progression géométrique de raison 5 

Si l’on veut mesurer la constante G par application des relations (3) et (13), 
la méthode la plus précise me paraît être l’observation des élongations néga- 


tives d’où l’on déduira a. Il sera nécessaire d'introduire des corrections de 


trois ordres, provenant, les premiéres de l’amortissement du système oscillant 
dont le décrément logarithmique sera désigné par d, les deuxièmes, du fait 


Fr l'équation (2) n’est pas linéaire à cause de la présence de Pinconnue a — 


ans son second membre; les termes correctifs correspondants font intervenir 
la déviation statique a, et se calculent au moyen des équations (11). Les troi- 
sièmes correspondent à Pexistence des harmoniques supérieurs dans l’expres- 
sion (1) du couple C, par l'intermédiaire des coefficients c,. 

Pai obtenu, de manière trop laborieuse pour que le détail puisse en être 
donné ici, les termes nécessaires pour pousser la précision jusqu’au dix- 
millième en admettant une déviation statique a, de l’ordre du centième de 
radian et un décrément logarithmique inférieur à 0,1, ce qui-est très facile- 
ment réalisable. Il sera nécessaire, pour éviter les perturbations d’origine 
électrostatique, d'employer un fil métallique pour la suspension. Avec une 
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durée d’oscillation de l’ordre de 100 secondes, il ne sera pas nécessaire d'aller 
au delà de la quatrième ou cinquième élongation négative après la mise en 
rotation du système B pour obtenir une déviation importante supérieure au 

dixième de’ radian. L’élongation négative de rang k est donnée par : 


1—e-td 1% Ihr / af x 
Ark Er (1-3) —35 G-3) 


= f e 2 
E (Es)] 


L’entier k wira pas au delà de 5 et æg se déduira de cette équation par une 
série facile et rapidement convergente d’approximations successives. Le terme 


principal du second membre est krag. 
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